PROCEEDINGS 

OF SCIENCE 



o 

O 

(N 

Or 

c/3 

(N 

(N 
> 
00 
(N 
O 
On 
O 

\o 

O 



X 



Investigation of the overlap of excited bottomonium 
states with hybrid operators 



C. Ehmann , T. Burch and A. Schafer 

Institutfiir Theoretische Physik 
Universitat Regensburg 
D-93040 Regensburg, Germany. 



E-mail: christian . ehmann@physik . uni-regensburg . de 



We analyze the overlap of color-octet meson operators with the T and the rji, and their excited 
states, especially the first radial excitations. Our analysis is based on NRQCD and includes all 
terms up to order v 4 . We use a variety of source and sink operators as a basis for the variational 
method, which enables us to clearly separate the mass eigenstates and hence to extract the desired 
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1. Introduction 

We present a study of couplings to the excited states of bottomonium through the use of a 
variational method. There are two reasons for this endeavor: We first wish to correct the interpre- 
tation given in Ref. [|l|] of there being no sign of low-lying excitations (e.g., radial) in heavy hybrid 
correlators by showing that they are indeed there after all, but that the coupling is simply too small 
to be accessible by the methods used in Ref. [jj]]. The point is that interpolators of any given quan- 
tum numbers will project out all states with those quantum numbers, regardless of whether the 
interpolators are hybrid in nature (i.e., containing gluonic fields; see, e.g., Ref. Our second 
reason is to demonstrate the efficacy of the variational method for finding couplings of excited 
hadronic states. This has been recently studied for excited pions [f||. Here, we show that ratios of 
couplings to a given state may be determined with good precision, even for resonance states and 
highly non-trivial currents. 



2. Method 

In the heavy quark limit we can readily define what a hybrid is: A hybrid excitation is a state 
which has large contributions of configurations where the quark-antiquark pair form a color octet. 
Such a qq pair is accompanied by an arbitrary number of valence gluons to give a singlet. Since 
these states, where the qq pair is in a defined color representation, are not mass eigenstates of 
QCD, mixing between such configurations will occur. So we expect hybrid operators to have finite 
overlap with all excitations of a given state. We show this for the bottomonium system, where the 
mixing is actually suppressed compared to mesons containing lighter valence quarks. 

To extract masses and couplings of excited states we rely on the variational method developed 
by Michael [Q] and later refined by Liischer and Wolff [Q]. 

In order to obtain the correlators, we first have to calculate the fermion propagators. Since 
bottom quarks are considered, the framework of NRQCD is applicable. We include all terms up to 
0(v 4 ), where v is the velocity of a quark, according to the power counting in 

The propagation of the fermions is given by: 

Hy,t + a) = (l - aJ ^^-) n utit) (l - -aSHmM, M 
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The tildes denote improved versions of the corresponding derivatives. We use n = 2, which is more 
than sufficient in our case. B and E are the magnetic and electric fields created via the usual clover 
formulation. 



The last two terms of ( |2.3| ) are responsible for the configuration mixing mentioned above. 

We determine the quark mass for our simulation from finite momentum correlators for the T 
by tuning the kinetic mass extracted from the non-relativistic energy-momentum dependence to the 
experimental mass of the T. Since we want to investigate the rjb and the T, we use pseudoscalar and 
vector currents. For both we have a normal and a hybrid version. Table ||] gives an overview of the 
local operators we use. The P-wave states are only needed to set the scale. In order to assemble our 
basis with more linearly independent operators, we additionally smear the quark and the antiquark 
field independently with two different smearing levels. In total twelve different operators, each at 
the source and the sink, are available for constructing the cross correlator matrix. 
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Table 1: Overview of the used operators (9,. The hybrid versions of the P- waves are not needed. 



We note in passing that the hybrid operators we use are related to higher twist contributions to 
the second moment of parton distributions (see [0] for a related topic): 

yEijkFjky ~ \ffeijk[Dj,D k ]Y. (2.4) 

We do not explicitly calculate such a contribution (for bottomonium such a quantity is probably of 
little interest), we just point out the usually inherent difficulty behind what it is that we are after, 
especially for the excited states. 

We determine the correlator matrix 

Qj(t) = (di(t)dj(0)) (2.5) 

and solve the generalized eigenvalue problem: 

Qj(0V»/ = Ki(t,to)Qj(to)Wj, (2.6) 

where A„ is the eigenvalue corresponding to the eigenvector iff". 

Since due to fluctuations Cu is not exactly symmetric (although it is within errors), we sym- 
metrize it by hand in order to make the diagonalization procedure more stable. 

The eigenvalues are given by [Q, ||] 

K(t,to) = A n e- M » {t - t0 \l + O(e- m ^-^)}, (2.7) 

where M„ denotes the mass of the nth state and AM n the mass difference to the next state. So for 
large enough values of t, we have a single mass state in each channel. 
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The eigenvectors shed light on the couplings of the states to the operators in use. Again, this 
should hold for large values of t. That is easy to see by regarding the following relation, where the 
eigenvalues are reconstructed by applying the matrix on the eigenvectors: 

Cij(t)Y] (OilVjOj) (Oj\n) 

QKOvf {OkWfOt) {o k \ n y 

where (0/|n) stands for the overlap of the ith operator with the nth eigenstate. Consequently, we 
are able to make statements about the ratios of couplings of different operators to the same state. 



3. Results 

We are working on configurations provided by the MILC -collaboration Q|. They were gener- 
ated using improved staggered fermions and the Liischer-Weisz gauge. Table || shows the param- 
eters of the lattices used. For the lattice spacing there are two values given. The first one comes 
from the analysis of the spin-averaged IP- IS splitting, the second one is given by the MILC- 
collaboration. 





a- l [MeV] 


volume 


N f 


am sea 


ami, 


8.40 


2378/2279 


28 3 x 96 





OO 


1.7, 1.8 


7.09 


2097/2252 


28 3 x 96 


2+1 


0.0062/0.031 


1.7, 1.8 


6.76 


1495/1587 


20 3 x 64 


2+1 


0.01/0.05 


2.4, 2.5 



Table 2: The three different lattices we use in our simulations together with the corresponding parameters. 
The first value for the inverse lattice spacing comes from the IP- IS splitting, the second one is given in 

Let us start with a five dimensional basis containing the operators Nll(l), Nln(2), Nnn(3), 
Nww(4), Hll(5). Here, the capital letter denotes the type of the operator (N=normal, H=hybrid), 
the lower case letters characterize the smearing of the quark and the antiquark (l=local, n=narrow, 
w=wide) and the number in brackets is the index of ; . 

Figure |l] shows the effective masses of the five eigenvalues of the T in this basis for the dynam- 
ical lattice with /3 = 7.09 and am/, = 1.7. The ratios of the couplings of the two local operators to 
the ground state and to the first radial excitation are plotted in Figure ||. One can see that this ratio 
is about 1/90. This is somewhat suprising, since one would naively expect that radial excitations 
have a larger gluonic content than the ground state. About the same ratio is obtained also for all the 
other radial excitations. However, the scenario changes for the fifth state, which we identify with 
the "hybrid excitation." In the left plot of Figure || one can see that the absolute value of the ratio 
of the couplings is about one. 

An important check of our analysis is that the ratio of the amplitudes of the local operators 
should be independent of the number of extended operators contained in the basis. The right plot 
in Figure |3| shows that, in fact, this is the case. The ratio is shown in the five dimensional basis 
given above and in the three dimensional basis Nll(l), Nnn(2), Hll(3). The same independence 
is also observed for all excited states. Therefore, we can conclude that the local operators are 
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Figure 1: Effective masses of the five eigenvalues in the basis Nll(l), Nln(2), Nnn(3), Nww(4), Hll(5) for 
the dynamical lattice with /3 = 7.09 and amt, = 1.7. to is the normalization timeslice. 

"approximately" orthogonal to the smeared ones. A similar sanity check also works for the large 
coupling ratio found for the hybrid excitation since it does not only appear as the highest state when 
we add more hybrid operators. 

Since lattices with different spacings are accessible to us, we can search for possible scale 
dependencies. The ratios of the couplings for the ground and first excited state for two different 
lattice spacings are shown in Figure |J A clear dependence is visible. Given that the hybrid oper- 
ator (in the numerator of the ratio) is not strictly a local current, but rather extended over a 2 x 2 
clover, this seemingly strong scale dependence is not that surprising. Unfortunately, with only two 
different lattice spacings available (and the same extent of the hybrid operator in lattice units) and 
outstanding renormalization we cannot say anything more definitiv about this. 

Similar results show up for the r\b, where the ratio for the non-hybrid states has changed to 
about 1/30, which can be traced back to the fact that in the h operator all three components of 
the fi-field are included. 
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Figure 2: The ratio of the couplings of the two local operators to the ground state is shown in the left plot, 
the one to the first excitation in the right plot. In both cases we use the five dimensional basis of Fig. 1. 



o 

-0.2 
-0,4 
-0,6 
-0.8 

-1 
-1,2 
-1,4 
-1,6 
-1.8 

-2 



P=7.09, 28x96 , N,=2+l, am b =1.7, t Q =l 



o cV./cV. 



€> » © 



J i I i L 

2 4 6 

t/a 



I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 



□ C. W / C. .W . 3-dim basis 
- 1 j t j ij T j 



i I i I 
8 100 2 4 6 



- 0,0115 



0,011 



0,0105 



0,01 

8 10 12 14 16 18 20 
t/a 



Figure 3: The ratio of the couplings of the two local operators to the hybrid excitation in the new basis is 
shown in the left plot. The right plot shows the ratio for the ground state in the three and five dimensional 
basis (the datapoints for the latter one are shifted by 0.3 in t/a for the sake of clarity). 
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Figure 4: The ratio of the couplings of the two local operators to the ground state is shown in the left plot, 
the one for the first excitation on the right plot. Each plot is given for two different lattices with = 2101 
MeV (fine) and a -1 = 1497 MeV (coarse). 
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